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MOMENTUM  FLUX  INCREASES  AND  COHERENT-STRUCTURE 
DYNAMICS  IN  A  SUBSONIC  AXISYMMETRIC  FREE  JET 

1.  INTRODUCTION 

Numerical  investigations  have  been  used  to  study  both  planar  and  axisymmetric  shear  layers.  A 
number  of  these  have  modeled  temporally  developing  mixing  layers  (Patnaik  et  al.  1976;  Riley  &: 
Metcalfe  1980;  Moin  et  al.  1985).  Others  have  modeled  spatially  developing  layers  (Ashurst  1979; 
Davis  &  Moore  1985;  Grinstein  et  al.  1986),  which  use  inflow  and  outflow  boundary  conditions  and 
therefore  correspond  more  closely  to  the  physical  situations. 

Only  in  very  few  cases  have  the  simulated  flows  been  compared  to  experimental  data  to  enhance 
our  basic  understanding  of  large-scale  organized  structures  and  turbulence  physics.  One  such  study 
involving  direct  numerical  simulations  of  the  plane  mixing  layer  using  spectral  methods  (Metcalfe 
et  al.  1987)  shows  good  agreement  with  experimental  data  on  educed  coherent  vorticity,  strain  rate 
and  production,  in  spite  of  the  clear  differences  between  the  simulated  and  measured  flows:  the 
simulations  addressed  a  spatially  periodic  temporally  developing  transitional  mixing  layer,  while 
the  laboratory  measurements  were  made  in  a  spatially  evolving,  temporally  stationary  mixing  layer 
of  a  much  higher  Reynolds  number. 

Recent  numerical  simulations  of  unforced,  subsonic,  compressible,  spatially  evolving  two- 
dimensional  shear  layers  have  shown  a  number  of  important  features  of  the  flow  field  in  the  transition 
region  (Grinstein  et  al.  1986,  1987).  The  calculated  distribution  of  merging  locations  (Grinstein  et 
al.  1987)  was  in  good  agreement  with  the  experimental  distributions  obtained  with  weakly  excited 
axisymmetric  (Zaman  &  Hussain  1977;  Kibens  1979)  shear  layers.  The  results  indicated  that  the 
presence  of  an  underlying  degree  of  organization  in  unforced  two-dimensional  shear  layers  is  a  result 
of  upstream  feedback  of  perturbations  caused  by  downstream  vortex  interactions. 

This  paper  represents  an  experimental- numerical  collaboration,  prompted  by  the  experimental 
result  that  the  momentum  flux  in  a  jet  was  found  to  increase  in  the  downstream  direction  ;  Clark 
1974;  Hussain  &  Clark  1977),  contrary  to  traditional  expectations.  The  result  and  interpretations 
presented  here  are  not  only  a  successful  result  of  the  collaboration  but  are  also  a  dear  evidence  of 
the  crucial  role  of  computer  simulations  in  answering  nuestinns  which  experiments  alone  can  not 
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easily  answer.  We  strongly  feel  that  such  collaboration  between  the  two  complementary  approaches 
is  necessary  for  understanding  the  physics  of  turbulence.  In  this  paper  we  investigate  the  validity 
of  the  traditional  hypothesis  of  isobaric  free  jet  flow.  It  is  usually  assumed  that  the  mean  static 
pressure  gradients  in  a  free  jet  are  everywhere  negligibly  small  compared  to  other  mean  forces  on  the 
fluid.  Thus  the  pressure  in  a  jet  is  considered  virtually  constant  and  equal  to  the  ambient  pressure 
(Abramovich  1963).  As  a  consequence  of  this  hypothesis,  the  mean  longitudinal  momentum  flux 
should  be  a  constant  at  each  streamwise  station.  Historically,  the  momentum  flux  in  a  free  jet  has 
been  considered  a  fundamental  invariant.  In  fact,  free  jet  data  which  failed  to  show  constant,  of 
momentum  flux  were  generally  rejected  as  unreliable  (e.g.,  Mons  et  al.  1971;  Kline  et  al.  1973). 

In  the  experimental  study  of  planar  jets  by  Hussain  k  Clark  (1977),  the  streamwise  momentum 
flux  was  computed  from  velocity  data  measured  with  hot  wires.  It  was  found  that  the  momentum 
flux  increases  in  the  streamwise  direction,  and  that  the  extent  of  the  increase  depends  on  the  initial 
conditions.  They  also  noted  that  pressure  can  drop  considerably  within  the  turbulent  regions  of 
a  free  jet  because  of  transverse  turbulent  fluctuations.  Although  they  measured  decreases  in  the 
static  pressure  that  were  consistent  with  the  observed  increases  in  the  momentum  flux,  they  could 
not  obtain  the  required  quantitative  balances  between  the  integrated  pressure  and  momentum  flux 
from  the  experimental  data.  This  was  attributed  to  large  uncertainties  in  both  velocity  and  pressure 
measurements  due  to  effects  of  high  turbulence  levels  in  the  jet,  especially  at  its  outer  edges,  on  both 
hot  wires  and  pitot-static- pressure  probes.  So,  while  the  streamwise  increase  in  the  momentum  flux 
in  a  free  jet  was  asserted  by  the  authors,  the  accuracy  of  the  measured  increase  in  the  momentum 
flux  and  of  the  measured  dependence  of  this  increase  on  the  initial  conditions  remained  unresolved. 

Here  we  use  time-dependent  spatially  evolving  simulations  of  shear  flows  to  study  the  transition 
region  of  a  subsonic,  compressible,  axisymmetric  free  jet  at  high  Reynolds  number.  In  addition  to 
having  the  advantage  of  providing  spatial  distributions  of  the  instantaneous  flow  field,  a  computer 
simulation  also  has  the  advantage  of  being  free  from  the  effect  of  probe  interference  and  the  highly 
restrictive  Taylor’s  hypothesis.  The  computations  have  also  the  advantage  of  being  free  from  the 
effects  of  high  turbulence  level  and  flow  reversal  on  hot  wire  measurements,  both  effects  character¬ 
istic  of  outer  regions  of  free  jet  flows.  The  calculations  show  that  negative  (i.e.,  countergradient  1 
production  occurs  in  regions  where  the  axis  of  an  inclined  structure  or  the  line  connecting  the  cen¬ 
ters  of  two  merging  structures  has  the  same  siope  as  the  mean  velocity  profile  (Hussain  1983).  This 
was  observed  in  laboratory  experiments  with  plane  mixing  iayers  subject  to  appropriate  forcing 
(Oster  k  Wygnanski  1982),  and  in  axisymmetric  jets  (Zaman  k  Hussain  1977,  1980)  in  the  region 
immediately  downstream  of  an  induced  localized  merging.  We  particularly  focus  on  the  question  of 
whether  the  momentum  flux  in  a  free  jet  flew  remains  constant.  The  '-alcufations  iuuicate  that  thr 
mean  static  pressure  drops  significantly  within  the  jet.  accounting  for  the  growth  in  the  streamwise 
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momentum  flux  and  supporting  the  experimental  results.  We  also  describe  the  spatial  features  of 
time-dependent  velocity,  vorticity  and  pressure  fieids  associated  with  coherent  structures  in  ine  jet 
and  discuss  their  dependence  on  controlled  excitation. 
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2.  THE  NUMERICAL  MODEL 


The  physical  system  we  simulate  is  a  compressible  axisymmetric  free  jet  of  air  emerging  into 
quiescent  ambient  air.  The  initially  laminar  jet  has  a  top  hat  axial  velocity  profile  at  the  exit  and 
uniform  standard  temperature  and  pressure  everywhere.  The  Mach  number  is  in  the  range  0.43- 
0.57,  the  free-stream  velocity,  U0,  is  1.5-2. 0  x  104  cm/s,  and  the  diameter,  D ,  is  1.4  cm.  The 
flow  system  is  shown  schematically  in  figure  1.  These  are  idealized  initial  conditions  since  most 
laboratory  jets  at  similar  Mach  numbers  are  not  initially  laminar.  The  jets  considered  here  are 
nearly  incompressible,  to  the  extent  that  the  turbulence  Mach  number  is  very  small.  At  the  same 
time,  the  Reynolds  number  and  the  Mach  number  for  the  fast-free-flow  regimes  considered  are  high 
enough  to  enable  their  optimal  study  with  an  inviscid,  explicit,  finite-difference  numerical  model 
such  as  discussed  below. 

The  numerical  model  used  to  perform  the  simulations  consists  of  the  two-dimensional,  time- 
dependent  inviscid  conservation  equations  for  mass,  momentum,  and  energy  density  for  an  ideal 
gas: 

!  =  -V.(,V),  U) 

=  -V  •  (pVV)  -  VP,  (2) 

BE 

—  =  -V.(EV)-V-PV.  (3) 

Here  E  =  P/( 7  -  1)  +  (1/2 )pV'i  is  the  total  energy  and  V,  P,  p,  and  7,  are  the  fluid  velocity, 
pressure,  mass  density,  and  the  ratio  of  specific  heats,  respectively. 

The  equations  are  solved  using  the  nonlinear,  high-order,  explicit,  compressible  finite-difference 
Flux-Corrected  Transport  (FCT)  algorithm  (Boris  and  Book  1976).  Through  a  two-step  predictor- 
corrector  scheme,  FCT  ensures  that  all  conserved  quantities  remain  monotonic  and  positive.  First, 
it  modifies  the  linear  properties  of  a  high-order  algorithm  by  adding  diffusion  during  convective 
transport  to  prevent  dispersive  ripples  from  arising.  The  added  diffusion  is  then  removed  in  an 
antidiffusion  phase  of  the  integration  cycle.  Hence,  these  calculations  maintain  the  high  order  of 
accuracy  without  resorting  to  artificial  viscosity  to  stabilize  the  algorithm. 

In  the  present  calculations  the  code  FAST2D  was  used  to  solve  equations  1-3.  The  code  uses 
direction  splitting  and  timestep  splitting  with  a  one-dimensional  fourth-order  phase-accurate  FCT 
algorithm.  The  code  FAST2D  has  been  extensively  tested  and  used  in  the  last  ten  years  (see  details 
and  references  in  Oran  Boris  1987).  Applications  to  subsonic  flows  include  calculations  of  vortex 
dynamics  and  asymmetric  mixing  in  shear  layers  (Grinstein  et  al.  1986.  1987).  The  simulations  of 
two-dimensional  compressible  shear-flows  have  reproduced  the  basic  large-scale  features  of  the  Mow 
observed  in  the  laboratory  experiments,  e.g.,  the  asymmetric  entrainment  (Grinstein  et.  al.  1986). 
the  distribution  of  merging  locations  (Grinstein  et  al.  1987),  and  the  spreading  rate  of  the  mixing 


layer  (shown  below  in  this  paper).  Much  work  has  been  done  to  test  boundary  conditions  and 
convergence  in  FAST2D  models.  Boris  et  al.  (1983.  1985)  described  tests  of  boundary  conditions  in 
both  subsonic  and  supersonic  flows,  and  Grinstein  et  al.  (1986)  described  boundary-condition  and 
resolution  tests  in  subsonic  and  choked  shear  hows.  In  addition,  Kailasanath  et  ai.  ( 1987)  reported 
tests  in  the  case  of  confined  subsonic  shear  flows  in  axisymmetric  ramjets. 

No  explicit  subgrid  turbulence  model  was  included.  However,  there  is  an  effective  filter  for 
high  frequencies  which  is  part  of  the  nonlinearity  of  the  FCT  algorithm.  The  filter  does  not  affect 
the  large-scale  structures,  but  ensures  that  wavelengths  smaller  than  a  few  computational  cells  are 
numerically  diffused.  Thus,  provided  that  the  computational  grid  is  chosen  fine  enough  to  resolve 
the  large-scale  features  of  the  flow,  the  residual  numerical  viscosity  of  the  algorithm  mimics  the 
behavior  of  physical  viscosity  at  high  Reynolds  number  by  dumping  small-scale  structures  on  the 
order  of  a  few  computational  cells.  There  are  essentially  no  viscous  effects  at  scales  greater  than 
four  or  five  cells. 

The  ability  of  the  nonlinear  filter  of  the  algorithm  to  simulate  the  effects  of  physical  viscosity 
was  verified,  in  particular,  by  a  study  of  the  laminar  spread  of  a  simulated  mixing  layer  by  Grinstein 
&  Guirguis  (1987)  which  will  be  reported  in  detail  elsewhere.  In  this  study  the  gridding  and  free 
stream  velocities  were  the  same  as  those  in  the  present  free  jet  flow  configuration,  but  zero  transverse 
velocity  was  enforced  in  the  calculations.  This  condition  (reasonable  in  the  high  Reynolds  number 
limit,  where  Vtrantv/Vstream  =  0(1/ Re))  was  enforced  in  order  to  ensure  that  the  spread  of  the 
shear  layer  was  strictly  laminar  and  due  only  to  the  residual  numerical  viscosity  of  the  algorithm. 
Then,  by  comparing  the  spread  of  the  velocity  profile  from  this  calculation  to  that  given  by  the 
full  numerical  solution  of  the  corresponding  incompressible  laminar  flow  equation  for  a  shear  layer 
(e.g.,  Schlichting  1968),  it  was  possible  to  obtain  a  bound  on  the  effective  numerical  viscosity  in 
the  simulations.  Using  the  bound  on  the  effective  numerical  viscosity,  and  based  on  the  jet  velocity 
and  diameter,  the  effective  Reynolds  number  Reo  for  the  large-scale  structures  in  the  present 
simulations  is  greater  than  as  1.5  x  10s. 

Inflow  and  outflow  boundary  conditions  are  used  to  represent  the  effects  of  the  region  outside 
of  the  computational  domain.  These  boundary  conditions  should  be  as  ‘‘transparent”  as  possible, 
which  means  that  unphysical  reflections  or  wave  absorption  at  the  boundaries  should  be  minimal. 
Imposing  such  boundary  conditions  is  difficult,  since  they  generally  require  representing  ‘he  flow 
field  outside  of  the  domain  with  a  very  restricted  amount  of  information  on  the  behavior  of  the 
solutions  at  the  boundaries.  Because  of  this  difficulty,  many  previous  simulations  have  used  periodic 
boundary  conditions  to  avoid  the  problem.  Such  temporally  >as  opposed  t0  sDatiaiivi  evolving 
calculations  can  be  thought  of  as  describing  the  time  evolution  of  the  flow  in  compact  regions. 


However,  laboratory  flows  evolve  both  in  space  and  time,  and  in  addition  to  the  limitations  of 
transformation  between  spatially  and  temporally  growing  disturbances  (  Gaster  1962:  Corcos  19801, 
the  dynamics  of  spatially  evolving  turbulent  shear  flows  cannot  be  captured  by  spatiailv  periodic 

simulations. 

The  calculations  presented  below  use  inflow  and  outflow  boundary  conditions  developed  and 
tested  for  multidimensional  FCT  calculations  (Boris  et  al.  1985;  Grinstein  et  al.  1986,  1987).  In 
order  to  compute  derivatives  at  the  boundaries  of  the  computational  domain,  we  define  the  behavior 
of  physical  variables  at  an  extra  row  of  guard  cells  surrounding  the  domain.  This  provides  a  way  to 
model  the  behavior  of  the  fluid  near  the  boundaries  in  terms  of  what  is  expected  of  the  flow  field 
outside  of  the  computational  domain.  In  the  calculations  presented  below,  the  inflow  boundary 
conditions  are: 

Pg  —  Pinflowi  (4d) 

=  ^inflow,  (46) 

Pg  =  Pu  (4  c) 

where  the  subscripts  g  and  1  refer  to  the  guard  cells  and  to  the  first  row  of  cells  inside  the 
computational  domain  (at  the  inflow),  respectively.  Thus  the  density  and  inflow  velocity  of  the  jet 
are  specified  and  the  energy  at  the  mflow  is  determined  from  a  zero-slope  condition  on  pressure. 
This  allows  the  pressure  at  the  inflow  to  vary  in  response  to  acoustic  waves  generated  by  events 
downstream.  An  important  result  of  allowing  this  feedback  is  that  the  instabilities  can  evolve 
naturally  in  the  calculation,  and  they  do  not  have  to  be  driven  continuously  to  be  reinitiated. 

The  boundary  conditions  at  the  outflow  guard  cells  define  the  density  and  velocity  through 
zeroth-  and  first-order  extrapolations  from  the  last  two  cells: 

Pg  —  Pn,  (5a) 

vg  =  '2vn  —  un_i .  (5b) 

These  boundary  conditions  are  consistent  with  fixing  the  inflow  density  and  velocity  of  the  jet.  The 
energy  at  the  outflow  guard  cells  is  defined  in  terms  of  the  density  and  velocity  given  bv  eqs.  i  5a- hi 
and  the  pressure  defined  by 

Pg  =  Pn  +  ~rr — rr-(P'imb  -  Pn),  (5c) 

\Ig  ij) 

where  Y  is  either  the  radial  or  axial  coordinate,  and  the  subscript  j  refers  to  the  location  of  the  lip 
of  the  nozzle.  Equation  5c  is  the  result  of  interpolating  between  the  pressure  value  at  the  boundary. 
Pn,  and  the  ambient  value  at  infinity,  Pam0.  It  ensures  that  the  pressure  equals  P.im0  at  an  infinite 
distance  from  the  nozzle.  The  siow  relaxation  of  the  pressure  towards  the  known  ambient  value  is 
necessary  because  the  pressure  at  the  inflow  is  not  specified,  but  calculated  from  other  variables. 
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By  giving  a  reference  pressure  value  we  avoid  secular  errors  in  the  calculations. 

A  more  realistic  transition  to  the  flow  region  of  interest  downstream  of  the  nozzle  (the  small 
dotted-box  region  in  figure  1)  is  obtained  by  surrounding  it  with  a  buffer  region  and  imposing  the 
inflow  and  outflow  boundary  conditions  on  the  outer  boundaries  of  the  computational  domain.  In 
addition,  a  portion  of  the  nozzle  is  included  inside  the  computational  domain.  The  effectiveness 
of  the  inflow  and  outflow  boundary  conditions  depends  on  the  dimensions  of  the  buffer  region  and 
depends  inversely  on  the  time  interval  over  which  the  simulations  are  carried  out.  The  effectiveness 
of  the  outflow  boundary  conditions  is  also  somewhat  dependent  on  the  nature  of  the  numerical 
dissipation  at  the  boundaries  which  will  determine  how  well  spurious  reflections  will  be  damped 
there.  The  practical  choice  of  the  computational  domain  and  the  gridding  was  motivated  by  these 
considerations. 

Figure  2  shows  a  schematic  diagram  of  a  typical  grid  used  in  the  calculations.  The  cells  are 
closely  spaced  in  the  radial  (r)  direction  across  the  shear  layer,  where  the  large  structures  form, 
and  they  become  farther  apart  as  the  distance  from  the  shear  layer  increases  (for  r  >  D).  The 
cell  separations  in  the  streamwise  direction  (x)  also  increase  in  size  moving  away  from  the  lip  of 
the  nozzle,  located  at  R  =  0.7  cm  and  x  =  x0  =  0.  The  cell  sizes  are  essentially  constant  in  the 
small  dotted-box  area  of  interest  in  figure  1.  In  general,  Courant  numbers  in  the  range  0.4  -  0.5 
were  used  in  the  calculations.  There  are  100  -  135  computational  cells  in  the  cross-stream  (radial) 
direction  and  330  -  410  cells  in  the  streamwise  (axial)  direction.  The  mesh  spacings  vary  in  the 
ranges  0.05  <  Ax  <  0.52  cm  and  0.02  <  A r  <  0.67  cm  for  the  coarser  grids,  and  in  ‘he  ranges 
0.03  <  Ax  <  0.29  cm  and  0.G1  <  A r  <  0.67  cm  for  the  finer  grids. 

In  the  case  of  initially  laminar  free  jet  flows  with  very  large  Reynolds  numbers,  shear  stresses 
and  dissipation  terms  are  significant  only  in  the  very  thin  layers  of  the  flow  on  the  surface  of  the 
nozzle  and  across  the  shear  layer.  The  nozzle  boundary  layers  are  important  in  determining  the 
flow  near  the  edge  of  the  nozzle  and  the  initial  thickness  of  the  shear  layer.  The  latter  thickness 
and  the  jet  diameter  provide  the  two  lengthscaies  for  describing  the  circulr*  jet  flow  in  the  initial 
region;  sufficiently  farther  downstream,  the  diameter  becomes  the  only  lengthscale.  In  the  free  shear 
flow  sufficiently  downstream  of  the  nozzle,  the  large-scale  coherent  structures  should  be  essentially 
independent  of  Reynolds  number. 

The  boundary  layers  in  the  nozzle  have  not  been  explicitly  included  in  our  calculations.  To 
determine  an  effective  initial  momentum  thickness  90  of  the  shear  layer  we  note  that  the  Strouhal 
number,  Sf*„  =  }90jU0 ,  associated  with  the  natural  instability  frequency  of  the  shear  layer,  lies 
in  the  range  0.0125  -  0.0155  (Husain  &  Hussain  1983).  Based  on  the  spectra  of  axial  velocity 
fluctuations  (see  section  3.1  below)  we  can  obtain  a  characteristic  value  =  0.014  in  this  range, 
if  the  initial  momentum  thickness  of  the  shear  layer  90  is  chosen  to  be  effectively  equal  t0  the 
thickness  of  the  rim  of  the  nozzle  in  the  simulations,  that  is.  one  cell  length  across  the  shear  layer. 


For  easy  reference  this  value  of  90  will  be  used  as  the  lengthscale  of  the  shear  layer.  In  order  to 
see  that  this  is  consistent,  we  can  use  a  typical  longitudinal  phase  velocity,  up  =  /A,  as  0.6 U0,  to 
estimate  90/\0  ss  St/ 0.6.  Thus,  for  St  in  the  range  0.0125  -  0.0155,  one  wavelength  (A0)  of  the 
natural  oscillation  of  the  shear  layer  downstream  of  the  nozzie  exit  is  resolved  with  39  -  48  ceils  in 
the  transverse  direction  and  16  -  19  in  the  streamwise  direction.  Therefore,  the  spatial  details  of 
the  large-scale  structures  and  their  spatial  evolution  are  well  resolved.  Moreover,  we  show  below 
that  there  is  good  agreement  between  the  calculated  and  experimental  results  when  this  choice  of 
&0  is  used  as  an  effective  lengthscale  of  the  shear  layer  in  the  simulations. 


3.  RESULTS  AND  DISCUSSION 


The  plan  of  this  section  is  as  follows.  In  Section  3.1  we  investigate  the  instantaneous  large-scale 
spatial  features  of  the  flow  with  special  emphasis  on  the  correspondence  between  the  vorticity  and 
pressure  fields.  Then,  in  Section  3.2  we  study  the  time-averaged  quantities  to  examine  the  effect  of 
the  initial  conditions  on  the  jet  flow.  Compressibility  effects  are  assessed  in  Section  3.3.  Coherent- 
structure  dynamics  are  studied  in  Section  3.4  in  order  to  understand  specific  observed  features 
of  the  flow.  The  streamwise  momentum  flux  of  the  jet  is  studied  in  Section  3.5  for  the  case  of 
unforced  jets  as  a  function  of  the  initial  shear-layer  thickness  and  Mach  number.  Finally,  the  effect 
of  controlled  excitation  at  the  jet  iiuet  is  brieflv  examined  in  Section  3.6.  The  sum  of  the  results 
presented  provides  a  detailed  picture  of  the  correspondence  between  the  increases  in  streamwise 
momentum  flux  and  the  formation  and  interaction  of  coherent  structures  in  the  axisymmetric  jet. 

3.1  Instantaneous  Flow  Visualization 

We  focus  here  on  results  for  the  case  of  a  Mach  0.57  jet  with  U0  =  2  x  104  cm/«.  and  0o  in 
the  range  0.01  -  0.02  cm.  Figure  3  shows  typical  instantaneous  contours  for  90  =  0.02  cm  (i.e.. 
I?/0O  =  70)  of  the  azimuthal  vorticity,  static  pressure,  and  velocities,  in  the  plane  of  symmetry  of 
the  jet.  Local  minima  of  the  static  pressure  occur  at  an  essentially  fixed  distance  xj  ss  ID ,  where 
the  vorticity  sheet  rolls  up.  The  vortices  move  downstream,  interact  with  each  other,  and  thereby 
spread  the  vorticity  until  the  central,  potential  core  of  the  jet  disappears  at  x.,  approximately 
between  4 D  and  5 D  (consistent  with  laboratory  observations).  The  observed  quasi-periodic  flow 
pattern  is  characteristic  of  axisymmetric  jets.  As  a  result  of  their  mutual  interaction,  the  vortex 
rings  merge  at  nearly  fixed  locations,  xt  as  2D  ,  X3  as  41?,  where  the  centers  of  the  vortex-ring  cores 
are  aligned  in  the  radial  direction  and  where  the  amplitudes  of  the  corresponding  subharmor'cs 
saturate  (Zaman  k  Hussain  1977;  Ho  k  Huang  1982 j.  The  first  radial  displacement  of  the  vortex 
rings  is  used  for  our  definition  of  the  merging  locations.  We  have  not  observed  the  vortex  rings  to 
roll  around  each  other  more  than  once  before  their  merging.  Thus,  this  gives  a  consistent  definition 
of  the  merging  locations  -  given  somewhat  upstream  of  the  locations  where  the  merging  structures 
become  indistinct  and  the  mergings  are  actually  compieted.  The  vortex  interactions  are  promoted 
by  small  displacements  of  the  newly  formed  structures  from  the  center  of  the  shear  layer  as  they 
are  convected  downstream.  These  displacements  are  likely  to  be  caused  directly  or  indirectly  bv 
pressure  and  velocity  fluctuations  originated  by  perturbations  downstream  at  vortex-ring  mergings 
and  at  the  end  of  the  potential  core.  The  vortex  rings  are  either  slowed  down  or  accelerated  in  their 
convection  downstream,  depending  on  whether  they  are  shifted  away  or  towards  the  centerline  of 
the  jet  by  the  radial  displacements.  When  the  displacements  are  properiy  phased  with  a  pair  of 
structures  sufficiently  dose  to  each  other,  the  structures  roil  around  each  other  by  mutual  induction, 
which  leads  in  turn  to  the  vortex  merging. 


Vortex-ring  merging  appears  to  be  the  mechanism  for  the  spread  of  the  shear  Saver  in  the 
present  simulations.  Spreading  due  to  viscous  diffusion  is  expected  to  be  much  less  important  in 
the  limit  of  high  Reynolds  number,  and  is  negligible  here  due  to  the  very  small  residual  numerical 
viscosity  of  our  inviscid  numerical  model.  Other  mechanisms  of  shear-layer  spreading  such  as 
continual  fractional,  partial  and  total  pairings  as  well  as  vortex  tearings  (Hussain  &  Clark  1981), 
were  not  observed  in  this  simulation.  In  this  paper  we  concentrate  on  the  transition  region  of  the  jet 
upstream  of  the  end  of  the  potential  core,  where  the  flow  dynamics  can  be  expected  to  be  captured 
realistically  by  the  simulation.  Three-dimensional  effects  not  considered  in  the  simulations  are 
definitely  important  near  the  end  of  the  potential  core,  so  that  any  discussions  of  two-dimensional 
structures  beyond  this  point  is  of  questionable  relevance  to  the  physical  situation.  However,  such 
two-dimensional  simulations  are  useful  for  comparison  with  future  three-dimensional  simulations 
in  order  to  assess  the  role  of  three-dimensionality  and  to  understand  the  differences  between  the 
flow  physics  in  the  different  dimensionalities. 

Figure  3b  shows  that  in  the  transition  region  of  the  jet,  the  pressure  drops  about  11%,  14%, 
and  22%  below  ambient  at  the  locations  of  vortex  roll-up  ,  and  vortex  merging  at  x\  and  i3, 
respectively.  Such  minima  are  associated  with  centers  of  passing  vortical  structures.  The  pressure 
rises  above  ambient  between  the  structures  because  of  the  occurrence  of  stagnation  points  there  in 
the  reference  frame  of  the  advected  structures.  Positive  pressure  changes  of  up  to  8%  are  calculated. 
The  pressure  distribution  in  space  is  a  direct  consequence  of  the  instantaneous  large-scale  coherent 
structures  in  the  flow.  The  correspondence  between  Instantaneous  pressure  and  vorricitv  field 
can  be  clearly  discerned  in  numerical  simulations  but  is  considerably  more  difficult  to  establish 
experimentally.  The  pressure  increases  are  discussed  in  detail  in  the  later  sections. 

The  spectra  of  the  axial  velocity  fluctuations  at  r  =  0,  the  center  of  the  jet.  and  r  =  D/2, 
the  initial  center  of  the  shear  layer,  are  shown  in  figure  4  at  selected  streamwise  locations.  The 
spectra  are  based  on  time  series  which  exclude  the  initial  flow  transient  and  have  a  time-duration 
of  28-56  vortex  roll-up  periods.  Just  downstream  of  the  nozzle  lip  (i  =  0.5D.D  and  r  =  D/2),  the 
spectra  shows  a  dominant  frequency  peak  at  ss  0.014  associated  with  the  vortex-sheet  roil-ep 
process.  The  spectral  amplitudes  corresponding  to  the  first  two  subharmonics  of  St/>„  grow  as  we 
move  downstream  and  attain  their  maximum  values  at  x  a  2D  and  x  s:  4D.  the  first  and  second 
merging  locations,  respectively.  A  similar  streamwise  evolution  of  the  spectrum  is  observed  at  the 
centerline  of  the  jet,  r  =  0.  As  noted  previously  fGrinstein  et  al.  1987),  the  presence  near  the 
exit  of  the  nozzle  of  the  spectral  peaks  corresponding  to  the  subharmonics  indicates  feedback  from 
the  downstream  events  on  to  the  inflowing  jet.  The  roil-up  frequency  ~  0.014  corresponds 
to  the  natural  instability  frequency  of  the  shear  layer  and  is  within  the  expected  range  0.0125  - 
0.0155  (Husain  and  Hussain  1983).  Note  that  experimental  values  of  Stq^  1  Zaman  A:  Hussain  1077) 


are  consistently  lower  than  the  theoretically  predicted  most  unstable  value  of  0.017,  the  difference 
attributed  to  feedback  effects  from  the  structures  downstream  (Hussain  et  al.  1986). 


Figure  5  shows  vorticitv  contours  in  a  plane  passing  through  the  jet  axis  for  the  simulation  with 
0o=O. 01  cm  (D/8 o  =  140),  one  half  the  value  of  80  in  the  case  of  figure  3.  Decreasing  the  shear-layer 
thickness  reduces  the  wavelength  of  the  initial  instability  and  hence  the  lengthscale  of  the  vortex 
rings,  which  allows  a  third  merging  to  occur  before  the  end  of  the  potential  core.  Defining  i,+1 
(i  >  1)  as  the  quasi-stable  location  of  the  ith-merging,  we  note  that  (x3  -  x3)  «  x2  for  both  jets. 
Moreover,  (x*  —  x3)  is  of  the  order  of  2  -  3  times  (x3  -  x2).  This  spacing  was  seen  in  previous 
simulations  of  two-dimensional  axisymmetric  jets  (Grinstein  et  al.  1987)  and  is  in  agreement  with 
experimental  results  obtained  for  weakly  excited  circular  jets  (Zaman  &  Hussain  1977;  Kibens 
1979).  It  is  also  in  good  agreement  with  the  predictions  of  a  feedback  formula  suggested  by  Laufer 
&  Monkewitz  (1980)  and  by  Ho  &  Huang  (1982),  that  describes  the  feedback  between  the  convection 
of  instabilities  downstream  and  the  upstream  propagation  of  acoustic  fluctuations.  The  upstream 
feedback  appears  to  be  responsible  for  the  underlying  degree  of  organization  in  the  axisymmetric 
free  shear  layer  (see  also  Hussain  &  Zaman  1978). 

3.2  Time-Averaged  Data 

Time-average  results  are  presented  here  for  the  unexcited  jet.  Although  these  results  do  not 
capture  the  dynamics  of  the  coherent  structures,  they  do  show  the  effect  of  the  localization  of  the 
roll-up  and  merging  processes  on  the  average  properties  of  the  flow.  Moreover,  these  results  will 
also  be  used  as  a  reference  to  evaluate  the  effect  of  compressibility  and  controlled  excitation  on  the 
jet. 

The  time  average  of  the  quantity  Q  is  calculated  based  on 


Q  =  Uni  ^  /  Q(t)dt , 

T-oo  T  J0 


where  the  integration  interval  is  chosen  in  such  a  way  that  excludes  the  initial  flow  transient,  and 
in  practice  we  use  a  finite  time  T  -  typically  in  the  range  of  32-72  vortex  roll-up  periods  -  which 
is  adequate  to  show  the  outstanding  statistical  features  of  the  results. 

In  figure  6  the  normalized  width  of  the  shear  layer,  B/80 ,  is  shown  as  a  function  of  x/80.  Here. 
B  =  (Rom  ~  #0.92)1  where  #*  is  defined  for  a  given  x  as  the  radial  location  at  which  7 I  =  f’  =  ki'}. 
The  figure  shows  the  approximately  linear  growth  of  the  mixing  layer  in  the  axial  direction  for  both 
calculated  and  experimental  results,  and  the  good  agreement  among  the  growth  rates  dB/dx.  The 
experimental  data  for  initially  turbulent  free  air  circular  jets  corresponds  to  the  case  with  lowest 
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inflow  fluctuation  level  studied  by  Hussain  &  Zedan  (1978),  and  to  the  experiments  of  Husain 
(1981).  In  the  simulations,  the  growth  rate  of  the  mixing  layers  decreases  immediately  downstream 
of  the  locations  of  quasi-stable  vortex  roll-up  and  vortex  merging  indicated  by  the  arrows  in  the 
figure.  The  relative  locations  of  the  mergings  depend  on  the  wavelength  of  the  most  amplified 
mode.  For  90/D  «  1,  as  we  have  here,  the  properties  of  the  axisymmetric  shear  layer  in  the 
near  field  of  the  jet  should  aproach  those  of  the  plane  mixing  layer.  Following  the  predictions  of 
linear  stability  theory  (Michalke,  1965),  this  wavelength  should  then  scale  with  90  for  a  given  jet 
velocity.  In  agreement  with  this  theory  figure  6  shows  that  the  distances  from  the  nozzle  edge  to 
the  vortex-ring  roll-up  and  merging  locations  scale  with  90.  Figure  7  shows  similar  variations  of 
the  axial  RMS  velocity  fluctuations  u' /U0  along  the  centerline  of  the  jet  as  a  function  of  the  ratio 
D/9oi  indicating  also  that  the  fluctuation  level  at  each  streamwise  location  x/90  roughly  scales 
with  90/D. 

Figures  8a-e  show  contours  of  the  time-averaged  streamwise  velocity  U /f/0,  Reynolds  stress 
puv/iPambUZ),  RMS  axial  velocity  fluctuations  u'/U0  =  (u2)*/U„,  static  pressure  (P~Pamb)/Pimb, 
and  RMS  cross-stream  velocity  fluctuations  v' /U0  =  (v2)% /U0,  where  u  =  u  —  U  and  v  =  v  -  V 
are  the  instantaneous  velocity  fluctuations.  Figure  8a  again  shows  the  approximately  linear  spread 
of  the  shear  layer,  with  the  maximum  spread  rate  of  the  flow  occurring  near  the  locations  of 
the  maxima  of  puv/(pamt,UZ)  (in  figure  8b)  and  those  of  u'/U0  at  the  center  of  the  shear  layer 
(in  figure  8c),  somewhat  upstream  of  the  corresponding  locations  of  vortex  roll-up  and  vortex 
merging  at  the  maxima  of  v'/U0.  The  spread  rate  decreases  downstream  of  the  third  maximum 
of  puv/(pambUo),  for  i  >  2D.  This  region  of  reduced  growth  rate  is  associated  with  a  region  of 
negative  Reynolds  stress  (discussed  further  below  in  Section  3.4). 

Figure  8d  shows  that  the  mean  static  pressure  varies  within  3%  of  ambient  pressure:  it  is 
generally  below  ambient  in  the  mixing  layers  and  above  ambient  in  the  potential  core  of  the  jet. 
The  locations  of  local  pressure  minima  coincide  with  the  maxima  of  v' /U0  in  figure  8e,  that  occur 
in  the  neighborhood  of  the  locations  of  vortex  roll-up  and  vortex  merging.  This  indicates  the 
crucial  role  of  the  large  scale  features  of  the  coherent  structures  whose  imprint  is  manifest  in  both 
the  cross-stream  velocity  fluctuations  increases  and  the  static  pressure  decreases  within  the  shear 
layer.  In  addition,  we  note  that  the  pressure  maximum  just  downstream  of  the  second  merging 
location  (also  discussed  below  in  3.4)  appears  clearly  associated  with  a  region  of  larger  u'  JU0  near 
the  centerline  of  the  jet  at  i  ss  2.3 D  in  figure  8c,  and  at  x  «  32O0o  in  figure  7.  in  the  neighborhood 
of  the  second  merging  location. 

3.3  Effects  of  Compressibility 

In  order  to  assess  the  effects  of  compressibility  on  the  flow,  we  have  studied  the  case  of  a  lower 
Mach  number  (A/=0.43)  unforced  free  jet  with  the  same  diameter  and  ratio  D/90=  140.  and  with 
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U0  =  1.5  X  104  cm/s.  Figure  9  shows  a  typical  instantaneous  flow  visualization  for  this  case.  For 
a  given  effective  initial  thickness  of  the  shear  layer,  the  effect  of  decreasing  the  jet  velocity  is  to 
increase  the  smallest  characteristic  lengthscale  of  the  jet,  namely  the  lengthscale  of  the  structures 
that  roll  up  with  the  natural  shear-layer  instability  mode.  This  is  expected  because  the  size  and 
spacing  of  the  rolled  up  structures  scale  on  90  and  UQ\  (in  a  given  experimental  facility  U0  and  90 
are  uniquely  related). 

Because  of  the  larger  structures  involved  in  the  flow,  the  spread  rate  of  the  shear  layer  due 
to  vortex-ring  roll-up  and  vortex-ring  pairings  is  higher  for  the  slower  jet  (see  figure  10)  and  the 
potential  core  becomes  shorter.  In  contrast  to  the  higher  Mach  number  case  discussed  above,  for  the 
same  ratio  D/9a,  where  up  to  the  third  merging  occurred  before  reaching  the  end  of  the  potential 
core,  only  the  first  and  second  vortex-ring  pairings  now  occur  (cf.  figures  5  and  9).  Figure  11  shows 
contours  of  mean  static  pressure  for  this  jet.  The  pressure  changes  are  also  relatively  smaller  than 
in  the  M  =  0.57  case,  consistent  with  the  decreased  magnitude  of  the  vorticity  (which  scales  with 
U0),  when  M  =  0.47.  The  peak  values  of  negative  pressure  within  the  shear  layer  region  and  of 
positive  pressure  within  the  potential  core  drop  about  30%  and  50%,  respectively.  Because  the 
ratio  of  negative  to  positive  peak  pressure  is  larger  in  this  case,  the  pressure-drops  in  the  shear 
layer  now  become  relatively  more  important  than  the  increases  within  the  potential  core.  This  is 
consistent  with  the  faster  spread  of  the  mixing  layer  and  results  in  a  faster  increase  in  the  streamwise 
momentum  flux  for  lower  Mach  number  (discussed  below  in  Section  3.6). 

3.4  Coherent  Structures 

In  this  section  we  examine  the  detailed  coherent-structure  dynamics  of  the  flow  in  order  to 
understand  the  origin  of  the  region  of  negative  Reynolds  stress  puv/{pambU~)  in  figure  8b  and  of 
the  mean  pressure  maximum  in  figure  8c.  We  first  examine  a  sequence  of  contours  of  instantaneous 
vorticity  and  Reynolds  stress  puv/(pamt,U%)  in  figures  12  and  13.  The  sequence  shows  some  of 
the  details  of  two  fairly  complete  cycles  associated  with  the  first  and  second  mergings.  The  tilt 
of  a  vorticity  distribution  with  an  elliptically  shaped  cross-section  can  be  defined  in  terms  of  the 
alignment  of  its  major  oxis  with  the  direction  of  strain  rate.  Likewise,  the  orientation  of  the  line 
connecting  the  centers  of  the  vortices  involved  in  a  pairing  process  can  be  used  to  define  the  tilt 
of  the  associated  elliptically  shaped  vorticity  distribution.  The  downstream  (positive!  : ilt  of  the 
line  connecting  two  pairing  vortices  (i.e.,  slope  opposite  to  that  of  the  velocity  profile)  shows  that 
there  is  a  preferred  momentum  transport  from  the  jet  towards  the  shear  layer  which  will  have 
the  effect  of  increasing  the  turbulent  kinetic  energy.  When  the  tilt  is  negative,  i.e..  when  the  line 
connecting  vortex  centers  is  aligned  with  the  mean  velocity  gradient,  the  momentum  dux  is  away 
from  the  mixing  layer,  and  hence  the  turbulent  kinetic  energy  decreases  <  Browand  &  Ho  1983: 
Hussain  1983;  Ho  &  Huerre  1984;  Zaman  &  Hussain  1977).  Figure  13  indicates  that  positive 


Reynolds  stress  (solid-line  contours)  predominates  at  the  stage  of  the  pairing  in  which  the  tilt  is 
positive.  Past  the  stage  when  the  vortices  are  aligned  vertically,  v/hicli  defines  the  pairing  location 
(provided  that  the  vortex  rings  roll  around  each  other  only  once  before  merging),  the  tilt  becomes 
negative,  and  negative  Reynolds  stress  (dashed  contours)  increasingly  predominates  in  the  region 
of  the  merging.  These  properties  of  the  Reynolds  stress  during  the  pairing  process  were  previously 
noted  in  experimental  studies  of  the  mixing  layer  of  a  circular  jet  (Zaman  &  Hussain  1977)  and  in 
temporal  simulations  of  a  plane  mixing  layer  by  Riley  k  Metcalfe  (1980). 

The  turbulence  production,  II  as  - TivdU/dr ,  measures  average  energy  transfer  from  the  mean 
to  the  turbulent  flow  when  -uv  and  dU/dr  have  the  same  sign  (II  >  0).  Otherwise,  counter¬ 
gradient  momentum  transport  takes  place  and  the  energy  is  transferred  on  the  average  from  the 
turbulent  to  the  mean  flow.  In  general,  if  succesive  pairings  occur  in  the  shear  layer  as  we  move 
downstream,  the  positive  tilt  predominates  in  the  flow  and  the  growth  of  the  mixing  layer  ensures 
that  it  remains  actively  turbulent.  Thus  the  regions  of  negative  II  will  not  be  important  for  the 
average  flow,  and  will  not  be  observed  provided  that  the  growth  of  the  shear  layer,  and  hence  the 
positive  tilt,  are  sustained. 

Negative  production  has  been  observed  in  plane  mixing  layers  where  the  growth  of  the  shear 
layer  has  been  suppressed  by  appropriate  forcing  (Oster  k  Wygnanski  1982).  By  localizing  succe¬ 
sive  mergings  in  axisymmetric  jets,  it  has  also  been  observed  in  the  region  immediately  downstream 
of  the  pairing  location  (Zaman  k  Hussain  1977;  Hussain  k  Zaman  1980).  In  the  present  simula¬ 
tions  we  observe  negative  time-averaged  Reynolds  stress  in  figure  8b,  immediately  downstream 
of  a  vortex  merging  occurring  at  x  as  2. ID  in  the  case  of  the  9a  =  0.01  cm  jet,  and  further  on 
downstream  (at  i  ss  4 D)  for  the  jet  with  90  =  0.02  cm.  In  this  region,  the  tilt  is  predominantly 
negative,  and  figure  6  indicates  that  the  shear-layer  growth  slows  appreciably.  Upstream  of  that 
location,  the  succesive  and  relatively  close  occurrence  of  the  initial  roll-up  and  the  first  and  second 
mergings  causes  a  steady  growth  of  the  shear  layer  with  predominately  positive  tilt,  and  hence  the 
preferred  positive  7>uv /{pambUl)  in  the  first  two  diameters  of  the  jet  in  figure  Sb.  However,  the 
third  merging  only  occurs  for  x  >  4 D  for  this  jet  and  this  explains  the  region  of  reduced  growth 
which,  in  conjunction  with  the  localization  of  the  second  merging,  determines  the  observed  negative 
time-averaged  Reynolds  stress.  Capturing  these  effects  associated  with  the  streamwise  growth  of 
the  shear  layer  through  succesive  vortex  mergings  is  possible  because  we  are  simulating  the  flow 
dynamics  for  a  shear  layer  evolving  both  in  space  and  time. 

Figure  14  compares  contours  of  instantaneous  vorticitv  from  the  present  simulations  with 
those  of  coherent  (  phase-averaged)  vorticitv  for  the  axisymmetric  jet  studied  by  Hussain  ,'v;  Zaman 
(1980).  The  experimental  data  corresponds  to  a  jet  excited  at  a  frequency  StD  =  0.85.  which 
induced  stable  periodic  mergings  at  x  =  i.7oD.  Stable  mergings  are  induced  in  the  simulations  by 


the  self-excitation  of  the  jet  via  upstream  feedback.  The  figures  show  that  there  is  a  considerable 
streamwise  extent,  immediately  downstream  of  a  stable  merging  location,  where  further  mergings 
are  suppressed  and  the  shear-layer  growth  is  reduced.  The  corresponding  regions  of  negative  pro¬ 
duction  are  shown  in  figure  15.  The  higher  peak  values  in  figure  15a  reflect  the  higher  degree  of 
organization  in  the  shear  layer  due  to  the  excitation. 

The  origin  of  the  observed  pressure  maximum  on  the  centerline  can  be  investigated  by  studying 
the  detailed  evolution  of  the  instantaneous  vorticity  and  pressure  fields.  The  sequence  of  frames 
in  figure  16  shows  a  cycle  of  the  second  merging  occurring  in  this  particular  jet  (D/90= 70)  for  x 
between  32?  and  5 D.  As  expected  and  noted  earlier,  the  pressure  drops  below  ambient  pressure 
in  the  cores  of  the  vortex  rings,  and  rises  above  ambient  pressure  between  the  rings  and  within 
the  potential  core  of  the  jet.  In  a  reference  frame  advected  with  the  vortex  center,  the  pressure 
will  have  a  minimum  at  the  vortex  center.  Similarly,  starting  from  the  vortex  center  the  pressure 
will  increase  outwards,  becoming  maximum  at  the  saddles,  which  are  local  stagnation  points  in  the 
moving  frame  (see  figure  17).  The  interaction  between  vortex  rings  leads  to  the  first  and  second 
mergings,  at  x  »  22?  and  x  »  42?,  respectively.  In  particular,  the  instantaneous  difference  between 
local  and  ambient  pressure  attains  maximum  values  within  a  toroidal  region  just  downstream  of  the 
first  merging  location  (frame  a).  As  the  pairing  processes  continue,  the  radius  of  the  high-pressure 
ring  initially  at  x  *  22?  (in  frame  a)  progressively  decreases  as  it  convects  downstream,  focusing 
on  the  centerline  of  the  jet  (frame  6),  at  x  «  4.22?  for  the  jet  with  D/90  =  70  (at  r  as  32?  for 
Dl&0  =s  140),  just  upstream  of  the  streamwise  location  where  the  second  merging  is  completed.  As 
shown  below  in  figure  24a,  which  shows  contours  of  P  for  the  jet  with  D/90  -  70,  the  mean  pressure 
has  a  region  of  local  maxima  in  the  region  surrounding  this  virtual  focal  point  on  the  centerline. 

Maxima  of  mean  pressure  on  the  centerline  have  not  been  reported  in  laboratory  experiments 
with  circular  jets.  Figure  18  shows  results  for  the  case  of  a  circular  free  air  jet  with  U0  -  2.0  x 
103  cm/s  and  2?  =  7.6  cm  (Toyoda  1986).  The  mean  pressure  decreases  monotonically  along 
the  centerline  from  values  above  to  below  ambient  pressure.  This  behavior  was  also  reported 
for  planar  jets  (Miller  &  Comings  1957,;  Hussain  &  Clark  1977).  The  observed  decrease  of  the 
pressure  values  on  the  centerline  of  the  jet  with  decreasing  Mach  number  (cf.  figures  8  and  11) 
suggests  that  these  disparities  between  numerical  simulations  and  laboratory  experiments  may  in 
part  be  due  to  compressibility  effects  (the  jets  in  the  experiments  were  much  slower  and  virtually 
incompressible).  Furthermore,  we  note  that  azimuthal  effects  are  not  accounted  for  in  the  present 
simulations  and  become  more  important  as  we  approach  the  end  of  the  potential  core.  The  possible 
role  of  azimuthal  effects  in  further  reducing  the  pressure  values  within  the  potential  core  can  not 
be  assessed  at  present. 
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3.5  Tests  of  the  Isobaric  Jet  Flow  Hypothesis 

The  mean  streamwise  momentum  flux  can  be  evaluated  from  the  calculated  primitive  variables. 
This  is  done  by  first  integrating  the  axial  component  of  the  momentum  balance  equation  (eq.  2) 
in  the  radial  and  azimuthal  directions  (with  the  boundary  condition  V  =  0  for  r  —  oo),  and  then 
time  averaging  the  resulting  expression.  The  result  is 


Jyoo _  yoo 

'  p[U2  +  tt2]rdr  +  /  Prdr  =  constant, 

o  Jo 


From  eq.  7  it  follows  that  if  the  jet  flow  is  isobaric,  the  streamwise  momentum  flux  is  constant. 
This  constancy  was  universally  assumed  for  many  years,  and  even  used  as  criterion  for  acceptance 
of  experimental  data.  However,  from  eq.  7  it  is  also  consistent  with  conservation  of  momentum  to 
have  an  increase  in  the  average  momentum  flux  as  long  as  it  is  balanced  by  appropriate  decreases  in 
the  mean  static  pressure  integral.  Measurements  of  such  pressure  decreases  were  reported  by  Miller 
&  Comings  (1957),  and  later  on  by  Bradbury  (1965),  Sami  et  al.  (1967),  and  Maestrello  k  McDaid 
(1971),  although  corresponding  increases  in  the  streamwise  momentum  flux  were  not  observed  or 
reported.  Variations  of  the  streamwise  momentum  flux  were  first  reported,  in  connection  with 
planar- jet  experiments,  by  Hussain  k  Clark  (1977),  who  found  an  increase  in  momentum  flux 
and  consistent  decreases  in  the  static  pressure  integral.  On  the  basis  of  the  momentum  equation 
integrated  in  the  cross-stream  direction,  they  argued  that  negative  changes  of  the  mean  pressure 
within  the  turbulent  region  are  supported  by  the  transverse  turbulent  fluctuations,  that  is,  ~P{y)  - 
P(oo)  =  ~pv2,  where  y  is  the  cross-stream  coordinate  in  the  planar- jet  geometry,  and  P(oo)  ~ 
P<xmb-  This  relation  between  negative  pressure  changes  and  velocity  fluctuations  is  apparent  in 
figure  6.  Hussain  k  Clark  showed  that  momentum  flux  increases  were  accompanied  by  decreases 
in  the  integrated  pressure.  However,  due  to  the  unquantified  effects  of  high  turbulence  levels  on 
both  velocity  and  static  pressure  measured  in  their  experiments,  the  authors  could  not  check  the 
momentum  balance  quantitatively.  The  experiments  also  showed  that  the  extent  of  the  momentum 
flux  increase  depends  on  the  initial  conditions. 

There  is  a  second  inherent  constraint  on  the  measurements.  Because  the  hot-wire  error  due 
to  effects  of  extremely  high  turbulence  levels,  including  flow  reversal,  is  unavoidable  at  the  outer 
edges  of  the  jet,  the  axisymmetric  cross-section  tends  to  make  the  area  integral  of  momentum  ( hence 
momentum  flux)  diverge  unless  the  integration  is  truncated  or  the  velocity  profile  constrained  to 
satisfy  V  =  0  for  r  —  oo  in  an  appropriate  way.  Numerical  simulation  is  a  useful  tool  to  investigate 
this  matter.  Within  the  accuracy  of  the  simulation,  the  integrals  of  mean  and  turbulent  velocities 
and  static  pressure  are  unaffected  by  large  turbulence  levels  (unlike  experiments)  and  the  integrals 
can  be  computed  much  more  precisely. 
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We  define  the  streamwise  component  of  the  total  average  momentum  flux  .VfT(x),  and  the 
pressure  integral  V(x),  as 


Mt(x) 


p{r,x,t){U{r,x)2  +  u(r,x,t)2)rdr, 


V(x)  =  f  P(r,x,t)rdr. 
Jo 


and  rewrite  eq.  7  in  terms  of  these  quantities,  using  the  values  at  the  nozzle  exit  as  reference, 

Mt(x)-Mt(x<>)  V{x)-V{x0) 

Mt(x0)  +  Mt(xo)  ‘  J 

Figure  19  compares  the  two  terms  of  eq.  10  for  the  unexcited  jet.  The  momentum  flux  increases 
up  to  9%  relative  to  its  value  at  the  nozzle  exit  and  exactly  balances  the  decrease  in  the  pressure 
integral.  This  balance  is  evidence  of  the  accuracy  of  the  numerical  simulation  of  the  velocity  and 
pressure  fields.  Comparing  figures  19  and  3,  we  see  that  fast  changes  in  Mt{x)  occur  in  the 
neighborhood  of  x3  as  D,x 2  as  2D ,  and  Z3  a  4 D,  locations  at  which  vortex  rings  are  relatively 
strong  and  well  defined  and  the  static  pressure  has  large-scale  local  minima.  Note  that  in  the  case 
of  a  less  organized  flow  field  where  the  pairing  events  are  not  localized  the  variations  of  Mr(x) 
can  be  expected  to  be  monotonic  as  in  the  data  of  Hussain  &  Clark  (1977).  More  disorganized 
unexcited  jet  flows  would  result  if  azimuthal  effects  were  included  or  if  the  run  times  (and  therefore 
the  averaging  times)  were  much  longer. 

Figure  20  shows  the  sensitivity  of  the  ratio  Mt{x)/ Mt(x0)  to  changes  in  the  cutoff  value  rmax, 
used  in  the  numerical  evaluation  of  the  integrals  in  eqs.  8  and  9.  The  growth  of  the  momentum 
flux  along  the  first  five  diameters  is  essentially  independent  of  rmax  for  rmox>  3 D.  The  cutoff 
values  used  for  the  radial  integrations  were  chosen  in  this  range.  The  typical  averaging  times  used 
to  evaluate  Mt{x)/Mt{x0 )  varied  in  the  range  of  56-72  roll-up  periods.  With  these  averaging 
times,  the  the  peak  values  of  ;Wt(j)/,Wt(x0)  at  x3,x2,  and  X3  were  estimated  to  be  converged 
within  2  -  6%,  depending  on  the  streamwise  location.  The  convergence  becomes  slower  as  we 
move  downstream,  e.g.,  at  x2  and.  x3,  since  fewer  merging  realizations  than  roll-up  realizations  are 
involved  in  the  averaging  process  for  a  given  averaging  time. 

Figure  21  compares  the  calculated  streamwise  momentum  flux  for  the  unexcited  flow  with 
the  available  experimental  data.  For  smaller  90,  more  vorticitv  is  concentrated  in  the  shear  layer, 
where  there  are  correspondingly  larger  decreases  in  the  mean  static  pressure.  According  to  eq.  9, 
the  contributions  to  V{x)  from  pressure  decreases  in  the  shear  layer  are  weighted  more  (through 
the  factor  r)  than  the  increases  within  the  potential  core  of  the  jet.  This  explains  the  larger  values 
of  Mt(x)/Mt(x0)  (up  to  for  33%  larger)  in  the  case  with  smaller  initial  thickness  and  lor 
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Ua  =  2  X  lO4  cm/s).  Moreover,  for  lower  Mach  number,  the  shear  layer  spreads  faster  (see  Section 
3.3)  and  the  momentum  flux  also  increases  faster,  in  agreement  with  the  observed  trend  in  the  case 
of  the  plane- jet  experiments  reported  by  Hussain  k  Clark  ( 1977).  These  authors  also  reported  the 
results  for  the  circular  jet  based  on  the  experimental  data  of  Crow  &  Champagne  ( 1971 ),  which  do 
not  appear  to  be  consistent  with  the  results  of  the  present  simulations.  However,  we  note  that  in 
contrast  with  the  case  of  the  plane  jet,  the  accuracy  of  the  velocity  measurements  in  the  outer  edge 
of  an  axisymmetric  jet  is  far  more  crucial  because  as  r  -+  oc  the  outer  part  of  the  jet  becomes  much 
larger  in  area.  Consequently,  hot-wire  data  cannot  be  used  reliably  to  evaluate  momentum  flux 
because  it  suffers  from  large  inaccuracies  at  the  high  fluctuation  levels  typical  of  the  outer  edges 
of  jets.  Other  experimental  methods  such  as  laser- Doppler  anemometry  and  particle  displacement 
velocimetry  can  overcome  this  aspect,  even  though  they  have  other  inherent  inaccuracies. 

3.6  Effect  of  Excitation 

In  this  section  we  study  the  sensitivity  of  the  results  to  the  presence  of  controlled  planar 
excitation  at  the  inlet.  Excitation  or  forcing  of  shear  layers  has  been  extensively  used  in  the 
investigation  of  the  basic  mechanisms  of  shear  flows  (e.g.,  Crow  k  Champagne  1971;  Zaman  & 
Hussain  1977,  1980;  Ho  k  Huang  1982).  In  the  present  calculations,  the  jet  was  excited  at  the 
inflow  by  imposing  a  sinusoidal  (coherent)  perturbation  of  fixed  frequency  f\  and  relative  amplitude 
ai  on  the  free  stream  inflow  velocity  U0,  i.e., 

U0  —  Ug  [1  +  o1sin(2x/1r)].  (12), 

which  provides  a  controllable  planar  excitation  of  the  jet  stream  at  the  nozzle  exit.  This  imposed 
perturbation  simulates  the  forcing  introduced  acoustically  in  laboratory  experiments  by  speakers 
located  in  settling  chambers  upstream  of  the  nozzle  exit  (e.g.,  Zaman  k  Hussain  1980).  In  addition, 
the  type  of  fluctuations  that  typically  occur  in  experiments,  such  as  free-stream  turbulence  and 
boundary  layer  fluctuations  in  the  nozzle,  were  also  simulated  by  introducing  random  perturbations 
of  the  streamwise  velocity  at  the  nozzle  exit  (Grinstein  et  al.  1987).  The  fluctuations  introduced 
included  lengthscales  ranging  from  the  initial  shear  layer  thickness  to  the  jet  diameter  (see  the 
Appendix).  The  mean  fluctuation  level  of  the  streamwise  velocity  of  the  unforced  jets  at  the 
nozzle  exit  was  less  than  0.5%.  In  the  absence  of  controlled  coherent  planar  excitation,  the  inflow 
fluctuations  affect  the  degree  of  localization  of  the  mergings  (Grinstein  et  al.  1987).  Differences  in 
the  time-averaged  results  due  to  inflow  fluctuations  with  this  level  were  very  small  and  essentially 
undistinguishable  from  the  unforced  case.  We  considered  low-level  coherent  i  1%  -  2%)  excitation 
of  the  jet  with  0o=O.O2  cm  at  the  nondimensionaiized  frequency  S(q  =  fxD(Ug= 0.98.  which 
corresponds  to  the  the  natural  shear-layer  mode  (associated  with  StH_  ssrt.014  for  this  jet.  i.e.. 


The  coherent  excitation  of  the  jet  controls  the  wavelength  of  the  initial  instability  by  pe¬ 
riodically  modulating  the  shear  layer  thickness  at  its  origin.  The  associated  periodic  vorticitv 
perturbation  propagates  in  the  shear  layer  which  acts  as  a  waveguide.  Excitation  controls  the 
lengthscale  of  the  structures  that  roll  up  in  the  near  field  of  the  jet  and  can  modify  the  turbulence 
intensity  level.  The  variations  of  the  streamwise  momentum  flux  are  closely  related  to  variations  of 
the  mean  static  pressure  and  longitudinal  turbulence  intensity,  and  thus  depend  on  the  particulars 
of  the  excitation.  As  noted  above,  there  is  a  correspondence  between  the  maxima  in  streamwise 
momentum  flux  and  the  minima  of  mean  static  pressure.  In  the  excited  flows,.  Mt{ x)  increases 
faster  near  the  nozzle  than  in  the  unexcited  case.  This  effect  is  associated  with  the  larger  and  more 
organized  roll-up  process  which  is  induced  nearer  to  the  nozzle  at  higher  frequencies. 

Figure  22,  shows  a  sequence  of  instantaneous  contours  of  vorticity  during  a  time  interval  equal 
to  four  excitation  periods.  Forcing  at  the  shear-layer  instability  frequency  enhances  the  organization 
of  the  flow.  The  high-frequency  excitation  forces  the  initial  vortex-ring  roll-up  to  occur  closer  to 
the  lip  of  the  nozzle.  Likewise,  the  merging  locations  have  also  moved  upstream:  the  first  merging 
occurs  at  x  «sl.3D,  and  1.9D,  and  the  second  merging  occurs  at  x  «3.1D,  as  compared  to  x  «2D 
and  x  «4D,  respectively  in  the  unexcited  case.  The  contours  of  Reynolds  stress  puv/(pami,Uj) 
(figure  23)  show  positive  peaks  upstream  of  the  locations  of  vortex  roll-up  and  merging  and,  as  in 
the  unexcited  case,  a  well  defined  region  of  negative  values  about  x  as  3.6D,  just  downstream  of  the 
second  merging  location.  The  values  at  the  first  two  positive  peaks  have  increased  relative  to  the 
unexcited  case  by  about  10%  and  60%,  respectively.  The  values  near  the  second  merging  location 
have  increased  by  about  90%  at  the  positive  peak  and  considerably  more  at  the  negative  peak.  The 
negative  peak  value  is  small  and  less  accurately  defined  for  the  unexcited  case. 

Comparison  of  the  mean  pressure  distributions  for  the  excited  and  unexcited  cases  shows 
qualitatively  the  same  patterns,  but  with  more  pronounced  local  minima  and  maxima  in  the  excited 
case.  The  mean  static  pressure  (figure  24)  again  shows  minima  in  the  regions  of  the  merging 
locations  and  a  relative  maximum  on  the  centerline  near  the  the  region  where  the  second  merging 
takes  place.  The  pressure  at  the  negative  peaks  vary  by  -2.4%  and  -3.0%  relative  to  ambient 
pressure,  as  compared  to  —1.8%  and  -2.0%,  respectively,  in  the  unexcited  case.  The  peak  values 
of  Reynolds  stress  and  mean  static  pressure  are  noticeably  larger  as  a  consequence  of  the  sharper 
organization  of  the  roll-up  and  pairing  process  caused  by  the  excitation.  This  is  also  reflected,  in 
the  behavior  of  the  longitudinal  turbulence  intensity  u ' /U0  along  the  centerline  shown  in  figure  25. 

The  growth  of  the  shear-layer  thickness  is  also  faster  (figure  25)  for  the  excited  jet.  although 
downstream  of  the  merging  locations  it  slows  down  appreciably.  The  faster  growth  of  the  mixing 
layer  in  the  excited  jet  is  reflected  in  the  growth  of  the  streamwise  momentum  flux  shown  in 
figure  25c.  The  streamwise  momentum  flux  has  pronounced  peaks  at  the  merging  locations.  The 
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narrower  peaks  are  a  consequence  of  the  better  localization  of  the  mergings.  The  momentum  flux 
now  increases  up  to  14%  of  its  value  at  the  nozzle  exit.  The  larger  increase  in  momentum  flux  due 
to  excitation  is  considerably  more  pronounced  in  the  region  of  vortex  roll-up  and  first  merging. 

We  thus  find  that  the  excitation  of  the  shear-layer  mode  organizes  the  jet,  enhancing  the  effect 
of  the  smaller  scales  present  in  the  unexcited  flow  (shown  in  figures  23-25).  Figure  25b  indicates 
that  excitation  widens  the  shear  layer  and,  as  a  consequence,  widens  the  jet.  Moreover,  the  figure 
also  shows  that  forcing  distorts  the  nearly  linear  growth  of  the  shear  layer  observed  in  the  unexcited 
case.  Excitation  of  the  shear-layer  instability  mode  enhances  the  organization  of  the  shear  layer.  In 
particular,  the  mergings  become  more  localized,  and  as  a  consequence  the  streamwise  momentum 
flux  is  increased  up  to  about  50%  above  the  results  in  the  unexcited  case.  The  Reynolds  stress 
increases  considerably,  particularly  near  the  second  merging  location,  where  the  observed  negative 
production  is  now  much  more  important. 
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4.  CONCLUSIONS 

We  have  presented  results  from  finite-difference  numerical  simulations  of  the  transition  region  of 
a  subsonic,  compressible  axisymmetric  jet  at  high  Reynolds  number.  We  have  studied  the  sensitivity 
of  the  results  to  changes  in  Mach  number,  initial  shear-layer  thickness,  and  the  presence  of  low- 
level  controlled  velocity  perturbations  at  the  jet  inlet.  We  have  chosen  the  excitation  frequency 
associated  with  one  of  the  characteristic  lengthscales  of  the  jet,  namely,  the  shear-layer  mode.  The 
following  conclusions  may  be  drawn  from  this  work. 

Self-excitation  of  the  two-dimensional  jet  due  to  upstream  feedback  induces  quasi-periodic 
mergings  at  nearly  fixed  streamwise  locations.  The  turbulence  production  is  predominantly  positive 
when  the  growth  of  the  shear  layer  is  sustained  by  the  succesive  relatively  close  vortex-ring  mergings. 
As  the  mergings  become  more  separated,  a  significant  decrease  in  the  growth  rate  of  the  shear 
layer  occurs  downstream  of  the  merging  location  with  consequent  negative  turbulent  production, 
indicating  that  there  is  energy  transfer  from  the  turbulent  to  the  mean  flow.  As  was  previously 
asserted  based  on  experimental  observations,  the  origin  of  this  effect  is  purely  kinematical. 

The  static  pressure  varies  significantly  within  a  free  jet  relative  to  that  of  the  surroundings. 
These  pressure  variations  are  associated  with  the  large-scale  coherent  structures  and  their  inter¬ 
actions.  The  static  pressure  decreases  within  the  shear  layers  and  increases  in  the  potential  core. 
Instantaneous  pressure  drops  of  the  order  of  20%  were  observed  at  the  locations  where  the  vortex 
rings  merge,  in  contrast  with  pressure  increases  of  order  8%  elsewhere  in  the  jet.  The  mean  static 
pressure  variations  ranged  within  2-3%. 

The  total  momentum  flux  in  the  transition  region  of  the  unexcited  jet  increases  above  the 
jet  exit  value  up  to  a  maximum  value  in  the  range  9  -  15%.  depending  on  Mach  number  and 
initial  conditions.  The  increases  are  exactly  balanced  by  decreases  in  the  mean  static  pressure 
integral,  which  is  consistent  with  momentum  balance  requirements  and  agrees  reasonably  well  with 
experimental  results.  The  streamwise  mean  momentum  flux  decreases  downstream  from  the  nozzle 
along  the  first  few  diameters  of  the  jet.  This  is  due  to  the  pressure  variations  within  the  potential 
core  of  the  jet,  and  particularly  to  the  occurrence  of  a  pressure  maximum  on  the  centeriine  near  the 
location  where  the  second  merging  is  completed.  As  a  consequence,  the  increases  in  total  momentum 
flux  in  the  transition  region  of  the  jet  are  due  to  the  turbulent  contribution.  This  is  not.  however, 
inconsistent  with  the  data  of  Hussain  Sz  Clark  i  1977)  within  the  first  four  slit-widths  of  streamwise 
extent,  where  the  turbulent  contribution  was  of  the  order  of  (50%. 

Because  this  numerical  simulation  resolves  sound  waves,  it  deals  with  finite  sound  speeds.  This 
requires  that  a  sound  wave  travel  no  more  than  one  computational  ceil  length  in  one  umestep.  and 
thus  small  timesteps  dictated  by  the  Courant  condition  are  used  in  the  calculations.  If  we  try  to 
simulate  the  relatively  slow  flows  investigated  in  the  experiments,  the  number  'if  nmestens  in  a 
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reasonable  calculation  is  considerably  large  and  the  computations  become  impossibly  expensive. 
Patnaik  et  al.  (1986)  have  developed  an  implicit  FCT  algorithm  that  can  bypass  the  Courant- 
condition  limitation  and  extend  the  use  of  high-order  monotonic  algorithms  to  allow  long  timesteps 
in  slow  flows. 

APPENDIX 

In  previous  numerical  simulations  of  axisymmetric  jets  (Grinstein  et  al.  1987),  we  studied  the 
effects  of  random  perturbations  in  the  inflowing  jet  stream.  In  experiments,  such  fluctuation  may 
be  caused  by  free-stream  turbulence,  primarily  due  to  unavoidable  acoustic  disturbances  upstream 
of  the  nozzle  exit  and  in  the  laboratory  and  boundary  layers  in  the  nozzle.  In  the  calculations,  the 
fluctuations  are  modeled  by  replacing  the  streamwise  velocity  U0  within  the  nozzle  at  the  inflow 
boundary  (for  r  <  R)  with  I70(l  +  p),  where  p  =  p(r,t)  is  the  perturbing  term.  The  term  p  =  p(r,t ) 
is  added  to  the  sinusoidal  term  in  eq.  12  when  coherent  excitation  is  also  present,  and  is  defined  by 
a  sine  Fourier  series  in  the  variable  R 

P(r,t)  =  jf  53  sin  +<*>"»)  ,  U.l) 

where  Sm(t)  is  a  time-dependent  amplitude  defined  by  (see  figure  26) 

2(t  tom)/^rrn  if  tom  S  t  <  tom  "t*  <5tm/2 

<*m(t)  =  2 (tom  -  t)/6tm  +  2,  if  tom  +  Stm/2  <  t  <  tom  +  Stm 

.  0,  otherwise. 

Each  perturbation  term  in  eq.  A.l  varies  between  zero  and  a  maximum  value  of  no  more  than  a 
fraction  p£,/A/  of  the  inflow  jet  velocity.  The  duration  Stm  (20-50  timesteps  %  0.05 r0  -  9.13  r0, 
with  r0  =s  ^o/(5fs„Cfo))i  the  maximum  amplitude  p°m,  and  the  phases  <pm(0  <  om  <  2t)  of  the 
terms  in  the  Fourier  series  are  randomly  generated  numbers.  Starting  at  t  =  tom  =  0.  a  set  f^, 
Pmi  <^mi»  is  generated  for  each  mode.  Each  set  of  random  numbers  is  fixed  until  t  =  tom  +  8tm,  at 
which  time  we  reset  tom  —  tom  +  Stm  and  update  the  set.  The  number  of  terms.  N.  is  chosen  to 
ensure  that  the  perturbation  includes  wavelengths  ranging  from  the  initial  shear  layer  thickness 
to  the  diameter  D  of  the  jet,  thus  including  the  interval  of  important  length  scales  of  the  system. 
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Figure  1 

Flow  configuration  for  the  simulated  axisymmetric  free  jet. 
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Figure  2 

Typical  grid  used  in  the  simulations. 
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Figure  1 3 

Instantaneous  Reynolds  Stress,  pi  u  -  H)(v  -  ~:)/{pambU~),  for  the  case  \l 
140:  A  =  0.00815..  full  line:  >  0.  dashed  line:  <  0. 
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Figure  14 

a)  Coherent  (phase- averaged)  azimuthal  vorticity,  normalized  with  StpUc/D,  Hussain  (1986); 

b)  Instantaneous  azimuthal  vorticity,  normalized  with  SteoU0/9OJ  present  simulations.  Stqa  = 
0.014,  M  =  0.57,  D/0O  =  140. 
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Figure  21 

Streamwise  momentum  flux  increase.  —  :  U0  -  2.0  x  104  cm/s  (M  =  0.57 ),D/90  =  140:  •  : 
U0  -  2.0  X  104  cm/s  ( M  =  Q.o7),D/d0  =  70:  O-:  U0  -  1.5  x  104  cm/s  (A/  =  O.43),Z>/0o  = 
140;  A:  planar  jet,  f70  =  0.8  x  103  cm/s,  Hussain  and  Clark  (1978);  ▲;  planar  jet.  U0  = 
1.5  x  103  cm/s,  Hussain  and  Clark  (1978);  a:  circular  jet.  (J0  =  3.8  x  103  cm/s.  Crow  and 
Champagne!  1971 ). 
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Figure  22 

Instantaneous  azimuthal  vorticity,  Q/f0,  A  =  1.3,  shear-layer  instability  mode.  M  = 
0.57  ,D/90  =  70. 
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